The axis of evil 
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We examine previous claims for a preferred axis at (b, I) » (60,-100) in the cosmic radiation 
anisotropy, by generalizing the concept of multipole planarity to any shape preference (a concept we 
define mathematically) . Contrary to earlier claims, we find that the amount of power concentrated 
in planar modes for £ — 2, 3 is not inconsistent with isotropy and Gaussianity. The multipoles' 
alignment, however, is indeed anomalous, and extends up to £ = 5 rejecting statistical isotropy with 
a probability in excess of 99.9%. There is also an uncanny correlation of azimuthal phases between 
I = 3 and I = 5. We are unable to blame these effects on foreground contamination or large-scale 
systematic errors. We show how this reappraisal may be crucial in identifying the theoretical model 
behind the anomaly. 



The homogeneity and isotropy of the Universe - also 
known as the Copernican principle - is a major pos- 
tulate of modern cosmology. Obviously this assump- 
tion does not imply exact homogeneity and isotropy, but 
merely that the observed cosmological inhomogeneities 
are random fluctuations around a uniform background, 
extracted from a homogeneous and isotropic statistical 
ensemble. One may expect that the ever improving ob- 
servations of CMB fluctuations should lead to the great- 
est vindication of this principle. Yet, there have been a 
number of disturbing claims of evidence for a preferred 
direction in the Universe [JEHSSISHlIHIlS; 
making use of the state of the art WMAP first year re- 
sults These claims have potentially very damaging 
implications for the standard model of cosmology. 

It has been suggested that a preferred direction in 
CMB fluctuations may signal a non-trivial cosmic topol- 
ogy (e.g. pi 0, 0, 0), a matter currently far from 
settled. The preferred axis could also be the result of 
anisotropic expansion, possibly due to strings, walls or 
magnetic fields |15|. or even the result of an intrinsically 
inhomogeneous Universe ^(|. Such claims remain con- 
troversial; more mundanely the observed "axis of evil" 
could be the result of galactic foreground contamination 
or large scale unsubtracted systematics (see 0, 0, 0, 
|2C| for past examples). A closer inspection of the emer- 
gence of this preferred axis is at any rate imperative. 

In addressing the issue of statistical isotropy it is of 
paramount importance to work in harmonic space, a mat- 
ter that has not always been fully appreciated. Consider 
a full-sky map, ^f(r). It is well known that the space 
of such functions is not an irreducible representation of 
the rotation group, i.e. it contains invariant subspaces, 
the so-called harmonic multipoles. One should therefore 
perform an expansion into Spherical Harmonic functions 
(eigenfunctions of L 2 and L z , where L is the angular mo- 
mentum operator) 
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a systematic examination of statistical isotropy can now 
be carried out. 

The close relation between statistical anisotropy and 
non-Gaussianity (loosely the emergence of structured 
shapes) has been extensively discussed in the litera- 
ture 0, • For an isotropic Gaussian process each 
ai m is an independent Gaussian random variable char- 
acterised by variance Ci. In each realization the power 
Ci is randomly, uniformly distributed among all ms, and 
therefore there is no evident "shape" to the multipole. 
One finds conflict with this prediction if, for a given mul- 
tipole £, one identifies an axis for which an uncanny 
amount of power is concentrated in a given m-mode, once 
the z-axis is aligned with n. In general reorienting the 
axis destroys this feature, so the onset of m-preference is 
also a measure of statistical anisotropy. 

For example planar configurations (see 0,0; E3) have 
a z-axis' orientation where all the power is concentrated 
in the m = ±£ modes. They maximally break cylin- 
drical symmetry around the z-axis. On the contrary, 
cylindrically symmetric configurations ("Jupiters") have 
a z— axis orientation where all the power is in the m = 
mode, so that the sphere is divided into latitude bands, 
or zones, without any longitude variation in the temper- 
ature. Even in these extreme cases if we reorient the 
axis the power spreads to other m-modes and the m- 
preference is obscured. So the issue of m preference is 
closely linked with that of anisotropy. 

It was suggested in pjj that a possible statistic for 
anisotropy is 



max ^2 m2 1 



(2) 



and examine the problem multipole by multipole. Each 
multipole £ is an irreducible representation of SO (3) and 



the average of L? 7 maximized by an appropriate choice of 
z-axis; the outcome of this statistic is twofold: the value 
Ti and the axis n^. This statistic favours high ms and 
so it works well in searches for planarity (i.e. m — ±£). 
However the emergence of a preferred axis could come 
about from undue concentration of power in any m, a 
remark that is the main point of this letter. 
We consider the alternative statistic: 
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1 = 5 in galactic coordinates 
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1 = 5 in preferred frame 
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FIG. 1: The I = 5 multipole for a cleaned map in galactic co- 
ordinates (top) and aligned with (b, I) « (50, —91) (middle.) 
We have superposed the multipole vectors, and the chain link- 
ing them. For comparison we plotted the I — 3 multipole in 
its preferred frame. 

where Ceo = a^o 1 2 and Ci m = 2\a( rn \ 2 for m > (notice 
that 2 modes contribute for m^O). This statistic pro- 
vides three basic quantities: the direction ne, the "shape" 
mi, and the ratio rg of power absorbed by multipole rri£ in 
direction rig . The statistic is never ambiguous except for 
planar 1 = 1 configurations, which may be interpreted as 
m = 2 or m = 1 modes with axis n 2 rotated by 90°. As a 
convention we adopt the (I = 2, m = 2) rendition for this 
configuration. In addition to the shape, direction, and 



strength of a given multipole we may also define inter-£ 
quantities, such as the angle between preferred axes for 
adjacent multipoles. 

We have applied this method to a variety of maps: the 
cleaned and Wiener filtered maps of 1231: the internal 
linear combination maps (see also |36|): the power 
equalized maps of J24|. Taking cleaned and Wiener fil- 
tered maps, for £ = 2,3 our method promptly reproduces 
the results of Q, |2]J : we find planar multipoles aligned 
with (b, I) « (60, —100). For higher £ we find some sur- 
prising novelties. As an example, in Fig. ^ we plot £ = 5 
(top two panels). Galactic coordinates obscure the fact 
that the multipole is very approximately a pure m = 3 
mode aligned with (b,l) — (50,-91). This is revealed 
by our statistic and we have plotted the multipole when 
reoriented along its preferred axis. For comparison we 
plotted the £ — 3 mode aligned with its preferred axis. 

A complete application of this scheme to higher mul- 
tipoles is shown in Fig. [3 The multipole alignment 
found for £ = 2,. ..,5 does not extend to higher multi- 
poles (even though there is a rough alignment of multi- 
poles £ = 6, 11 around a strip on the b w 15 latitude). 
However it is highly inconsistent with isotropy. We per- 
formed Monte Carlo simulations for Gaussian maps with 
the best fit power spectrum, subject to the WMAP noise 
and beam. We then considered the average value of all 
possible angles between vectors for £ = 2, 5. This is 
of the order of 20° for the data. We find that only 5 out 
of 5000 simulations have an average value smaller than 
this. Thus we may reject isotropy on these scales at the 
99.9% confidence level. 

Other apparent anomalous features, however, are not 
as statistically significant as they appear. The ratios re, 
for instance, arc completely consistent with Gaussianity, 
as the top panel of FigEI shows. For example, for £ = 3 
we find that in the preferred frame 94% of the power is 
concentrated in the m = 3 mode. This was seen in pj 
as anomalous, at the 93% confidence level. Indeed this 
is true if we insist upon planarity. If on the contrary 
we allow for anisotropy due to any kind of m-preference, 
then the majority of the simulations actually find a m 
and n with a higher concentration of power in a single 
m; i.e. the observed rg is in fact low. 

Likewise, the apparent preference for low m multipoles 
(see Fig.|5Jl is not exceptional, once a more rigorous eval- 
uation of the feature is obtained from Monte Carlo sim- 
ulations (notice that the favoured ms are not uniformly 
distributed.) 

However we do find another unusual feature with some 
statistical significance. Phase correlations have been pro- 
posed as a discriminator for localized non-Gaussian fea- 
tures in the sky pfij . In our formalism we may define 
£ azimuthal phases by examining the phases of the ae m , 
for m > 0, as measured in the preferred frame. Most 
important will be the phase of ae m for the m that re- 
ceives most of the power in this frame. This may be used 
to fix the x and y axes of the preferred frame; because 
there is an overall sign ambiguity between different £ and 
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FIG. 2: A plot of m-preference results up to £ — 20. We 
seek the direction ne (galactic coordinates in the bottom two 
plots) where a given m (second plot from the top) receives the 
highest proportion of the power (with ratio ri plotted in the 
top panel). We have also plotted the Monte Carlo inferred 
variance for rg. We plotted results for cleaned (solid), Wiener 
filtered (dash), and internal linear combination (dotted) maps 
as described in text. 



m modes, phases separated by 7r should be identified. 
Inter-£ features may then be studied using the full set of 
Euler angles relating the preferred Cartesian frames for 
different multipoles. 

We find a very close alignment of the phases for the 
£ = 3 and £ = 5 modes. This should be blatant from 
Fig. ^ where we see that the longitudinal stripes of the 
two modes align very strongly. Comparison with Monte- 
Carlo simulations reveals that this feature is unusual to 
94.5%. Other phases don't correlate significantly. 

This concludes the presentation of our basic results, 
which may be summarised as follows. The planar shape 
of the low £ multipoles is not statistically significant, 
once we allow the simulations to seek an axis maximizing 
power in a single, but general m. However the alignment 
of the preferred axis is indeed very significant and extends 
up to £ — 5. An additional azimuthal phase correlation 
brings to 99.995% the confidence level for rejecting Gaus- 
sianity and statistical isotropy. 

We now examine the robustness of these results. 
Firstly they may be confirmed using an alternative tech- 
nique: the Maxwell multipolc vectors. These were first 
introduced by Maxwell |2^| over a century ago and 
have found a recent revival within the CMB commu- 
nity H H El E EJ HI H3. They encode the 2£ 
degrees of freedom contained in a multipole beside the 
power spectrum (see [2^, |3(], |^ and [3 for how these 
split into invariants and "eigenvectors"). Multipole vec- 



tors may potentially pick up m-preference since for pure 
Y£ m configurations £ — m multipole vectors align with 
the preferred direction n, with the remaining m spread- 
ing at equal angles in the orthogonal plane - we call this 
a "handle and disc" structure, and it provides a more 
visual approach to detecting m-preference. 

We consider these headless vectors in terms of pairs of 
points on the unit sphere. A handle is made through a 
clustering of points; a disc through a set of points trac- 
ing out an equator on the sphere. We look for these 
patterns by searching for the shortest routes between the 
points. The method does not provide a thorough search 
for m-preference, but can be used as a more visual tool 
to confirm and explore already noted features. In Fig. ^ 
we show how the method fares with £ — 3 where we find 
a single disk, and I = 5 where we find a nearly exactly 
planar disk of 3 vectors and a handle of 2. More gen- 
erally we confirm all the anomalous features described 
above. Where the two methods disagree there is no hint 
for m-preference. 

How dependent are these results on the chosen data- 
set, i.e. on the concrete rendition of the large-scale 
anisotropics after the galactic emissions have been re- 
moved? We find that most publicly available maps agree 
on the features outlined above and only differ signifi- 
cantly beyond £ = 12. The cleaned and Wiener filtered 
maps of [23j in fact agree almost exactly up to £ — 12 (see 
Fig. Some features below £ = 12, however, are not 
as robust when we consider other maps. The quadrupole 
features, for instance, may be easily erased in the rendi- 
tions of H|23 (see Fig. 0). The planarity of the £ = 6 
mode is present in some renditions (and is detected by 
our statistic) but not in others. Thus we should regard 
these features as more fickle. 

This leaves us with the obvious question: what could 
cause the strange directionality of the large angle CMB 
fluctuations? Most obviously they could be due to sys- 
tematic errors or galactic foregrounds. We subjected 
galactic templates 36] for synchrotron, free-free and dust 
emissions to this analysis (Fig. Unsurprisingly we 
found a strong non-Gaussian signal; however nothing cor- 
related with the direction or signal detected in the CMB 
maps. Specifically we found preference for m = modes 
for even £ and for m = 1 for odd £. This preference 
is invariably found at b — 90, i.e. with perfect align- 
ment with the galaxy axis. The synchrotron and dust 
maps have very similar morphologies; the free-free emis- 
sion is quite different and leads to more Gaussian m ra- 
tios. None of this correlates with the features claimed 
for the CMB, even when the templates are mixed with a 
Gaussian component. 

We also considered the simulated noise maps described 
in 37]. These include all known systematic effects. We 
found that none of the available 110 simulated maps dis- 
played any of the effects discussed in this paper. 

Should the observed preferred axis be real our remarks 
may be crucial in identifying the culprit theory. We found 
that although there is a preferred direction for low mul- 
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FIG. 3: m-preference results up to ^ = 20 for the V-band of 
galactic templates. We have included free-free (solid) syncro- 
ton (dot-dash) and dust (long-dash) maps; for reference we 
have reproduced results for the cleaned maps (dotted line). 

tipoles (and to some extent a preferred Cartesian axis) 
this does not pick a specific m. Hence we don't need a 
model favouring specific shapes (e.g. planarity); merely a 
model with a preferred axis. If the low £ fluctuations are 
due to the gravitational potential on the last scattering 
surface we can go further. The process may be described 



in terms of the potential Fourier modes $(k, 77^ ) ; each 
of these modes is reflected in the CMB in a pattern with 
exact m — shape, with n aligned with k. One must 
superpose several such modes to obtain a given £, specif- 
ically: 

a im = Aj2H^Vis)^MkA ms )Y; m (k) (4) 

Superposing a large number of modes leads to no pre- 
ferred direction or m preference. Should the number of 
modes be limited to a lattice, though, a preferred shape 
and/or preferred axis will emerge. For a cubic lattice the 
low £ may become superpositions of mainly three modes 
a em = A 1 5 m0 +A 2 Y* n (^ ,0)+A 3 Y* m (% , f ) and a solution 
in terms of Ai may now be found for a given observed 
m-preference. 

But it could also be that we live in a slab space, where 
there is a large number of modes in all but one direc- 
tion. This will erase any preference for a specific m, while 
keeping the preferred direction. The choice between the 
two possibilities hinges crucially on the phase correlation 
found for £ = 3 and £ = 5, with the implication that there 
may be a preferred frame (rather than just a preferred 
axis). We find this feature the most tantalising aspect of 
our analysis. 
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